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Abstract— Based on currently available kinematic data, we have searched for stars 
outside the Hipparcos list that either closely encountered in the past or will encounter in 
the future the Solar system within several parsecs. For the first time, we have identified 
two single stars, GJ 3379 (G 099-049) and GJ 3323 (LHS 1723), as candidate for a close 
encounter with the solar orbit. The star GJ 3379 could encounter the Sun more closely to 
a minimum distance d min = 1.32 ± 0.03 pc at time t min = —163 ± 3 thousand years. We 
have found two potential candidates for a close encounter that have only photometrical 
distances: the white dwarf SSSPM J1549-3544 without any data on its radial velocity 
and the L-dwarf SDSS J1416+1348. The probabilities of their penetration into the Oort 
cloud region are 0.09 (at a model radial velocity \V r \ = 50 km s _1 ) and 0.05, respectively. 

1 INTRODUCTION 

The Oort comet cloud (Oort 1950) is presumed to be located at the outer boundaries of 
the Solar system. It is highly likely that the cloud has a spherical shape and a radius 
of about 1 x 10 5 AU (0.485 pc). Close encounters of Galactic field stars with the Solar 
system play an important role in the dynamical evolution of the Oort cloud. In particular, 
such passages of stars can provoke the formation of comet showers that reach the region 
of the giant planets (Emel'yanenko et al. 2007; Leto et al. 2008; Rickman et al. 2008). 
Several researchers associate the traces of comet bombardments of the Earth with the 
impact of such showers (Wickramasinghe and Napier 2008). 

The question about the encounters of stars with the Sun within distances r < 2 — 
5 pc was considered by Matthews (1997) and Miillari and Orlov (1996) using various 
ground-based observations and by Garcia-Sanchez et al. (1999; 2001), Dybczyhski (2006) 
and Bobylev (2010) based on Hipparcos (1997) data in combination with stellar radial 
velocities. About 160 Hipparcos stars are known from the solar neighborhood 50 pc in 
radius that either encountered or will encounter the Solar system within r < 5 pc in a 
time interval of ±10 Myr. 

The goal of this study is to search for candidate stars closely encountering the Sun 
based on currently available kinematic data on stars that do not belong to the Hipparcos 
catalogue. Indeed, since the paper by Miillari and Orlov (1996), who analyzed all stars 
from the catalog by Gliese and JahreiB (1991), new observational data have appeared. 
In this paper, we solve the problem of statistical simulations by taking into account the 
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random errors in the input data in order to estimate the probability of a star penetrating 
into the Oort cloud region. 



2 THE METHODS 

2.1 Orbit construction 

We use a rectangular Galactic coordinate system with the axes directed away from the 
observer toward the Galaxy center (1=0°, b=0°, the X axis), in the direction of Galactic 
rotation (/ = 90°, b =°, the Y axis), and toward the North Pole (b = 90°, the Z axis). The 
corresponding space velocity components of an object U, V, W are also directed along the 
X, Y, Z axes. The epicyclic approximation (Lindblad, 1927) allows the stellar orbits to be 
constructed in a coordinate system rotating around the Galactic center. The equations 
are (Fuchs et al, 2006) 

X(t) = X(0) + ^ sin(^) + ^(1 - cos(k*)), (1) 

K Aid 

Y(t) = Y(0) + 2A(X(0) + - ^V(0) Bm(«t) + ^tf(0)(l - cos(«t)), 

Z(t) = ^® sin(z/t) + Z(0) cos(ut), 

where t is the time in Myr (we proceed from the ratio pc/Myr = 0.978 km s _1 ); A and 
B are the Oort constants; k = is the epicyclic frequency; Qq is the angular 

velocity of Galactic rotation for the local standard of rest, Qq = A — B; v = y/4nGpo is 
the vertical oscillation frequency, where G is the gravitational constant and p is the star 
density in the solar neighborhood. 

The parameters X(0), F(0), Z(0) and 17(0), V(0), W(0) in the system of equations (1) 
denote the current stellar positions and velocities. The displacement of the Sun from 
the Galactic plane is taken to be Z(0) = 17 pc (Joshi, 2007). We calculate the velocities 
U, V, W relative to the local standard of rest (LSR) with (U, V, W) LSR = (10.00, 5.25, 7.17) 
km s _1 (Dehnen and Binney, 1998). 

At present, the question about the specific values of the Sun's peculiar velocity relative 
to the local standard of rest (U, V, W) lsr is being actively debated (Francis and Anderson, 
2009; Binney, 2010). Arguments for increasing the velocity Vlsr from 5 km s~ x to 11 
km s _1 are adduced. Since, in our case, both the solar orbit and the stellar orbit are 
constructed with these values, the influence of the adopted Vlsr is virtually imperceptible. 
This is confirmed by good agreement between the encounter parameters of, for example, 
the star GJ 710 and other stars obtained with various values of (U, V, W)lsr by Bobylev 
(2010) and Garcia-Sanchez et al. (2001). 

We took p — 0-1 M Q /pc 3 (Holmberg and Flinn, 2004), which gives v = 74 km s _1 
kpc -1 . We used the following Oort constants found by Bobylev et al. (2008): A = 
15.5 ± 0.3 km s" 1 kpc" 1 and B = -12.2 ± 0.7 km s" 1 kpc" 1 ; k = 37 km s" 1 kpc" 1 
corresponds to these values. 
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Note that we neglect the gravitational interaction between the star and the Sun. 



2.2 Statistical Simulations 

In accordance with the method of Monte Carlo statistical simulations, we compute a set 
of orbits for each object by taking into account the random errors in the input data. For 
each star, we compute the encounter parameter between the stellar and solar orbits, d = 



y AX 2 (t) + AY 2 (t) + AZ 2 (t)). We characterize the time of the closest encounter by two 
numbers, d min and t min . The errors in the stellar parameters are assumed to be distributed 
normally with a dispersion a. The errors are added to the equatorial coordinates, proper 
motion components, parallax, and radial velocity of the star. We separately consider stars 
with spectroscopic distance estimates obtained with typical errors of 20-30%. In this case, 
the random errors are added to the distance during the simulations. 



3.1 Stars with Trigonometric Parallaxes 

First, let us consider the solar neighborhood about 10 pc in radius using the list of 100 
nearest stars from the Chilean RECONS site (http://www.chara.gsu.edu/RECONS/). 
This list reflects the results published before January 1, 2009. The trigonometric parallax 
of each star in the list was calculated as a weighted mean of the results of several (from 
one to four) observations. We are interested in ten RECONS stars that are not Hipparcos 
ones (see Table 1). Other RECONS stars are ten M dwarfs without any radial velocity 
estimates and the remaining stars are Hipparcos ones. 

The equatorial coordinates and proper motions of the stars in Table 1 were taken 
from Salim and Gould (2003) and Lepine and Shara (2005), where the proper motions 
were calculated by comparing the 2MASS stellar positions with those from the Palomar 
Observatory Sky Survey. 

We took the radial velocities from the works of various authors. Note that, according 
to the original first-class measurements by Nidever et al. (2002), two stars, GL 406 and 
GL 905, exhibit a high stability in an observing time interval of 1-2 yr. The internal 
measurement error is ±0.1 km s" 1 ; therefore, they were suggested as candidates for stan- 
dards to determine the radial velocities. Since these stars are M type dwarfs, the external 
error in V r is ±0.4 km s _1 . We took the radial velocity of the center of mass for the triple 
system GL 866 ABC from the plot in Delfosse et al. (1999). Note that, on the whole, 
the radial velocities for the stars listed in the table differ insignificantly from those in the 
catalog by Gliese and JahreiB (1991). Only the star GJ 473, for which the discrepancy 
in V r is two orders of magnitude, constitutes an exception. We will discuss this situation 
below. 

For each star, we constructed its orbit relative to the Sun in the time interval from -2 
Myr to ±2 Myr. Data on six stars with an encounter parameter d < 3 pc are presented 
in Table 2 and the trajectories of these stars are shown in Fig.l. 




3 RESULTS AND DISCUSSION 



3 



(1) The parameters of the GL 905 encounter with the solar orbit found here are in 
good agreement with the estimates by Miillari and Orlov (1996), d m i n = 0.95 ± 0.11 pc 
and t m in = 36.3 ± 1.4 thousand years; our values have considerably smaller random errors 
due to the currently available data being very accurate. 

Note that the encounter parameter for Proxima Cen is d min = 0.89 ±0.02 pc (Bobylev, 
2010) or d min = 0.95 ± 0.04 pc (Garcia-Sanchez et al., 2001), which was obtained with a 
slightly different radial velocity. As we see from Fig.l, GL 905 has a chance to be even 
slightly closer to the Sun than Proxima Cen in ?«37 thousand years. At present, we know 
only two stars with a smaller (than that of GL 905) encounter parameter: GJ 710=HIP 
89825 with d min = 0.31 ±0.17 pc (Bobylev, 2010) and HIP 85661 with d min = 0.94 ±0.71 
pc (Garcia-Sanchez et al., 2001). 

(2) We have identified the stars GJ 3379 (G 099-049) and GJ 3323 (LHS 1723) as 
candidates for a close encounter with the solar orbit for the first time. Their trigonometric 
parallaxes with a relative error e n /n < 1% were first determined only recently by Henry et 
al. (2006). Their radial velocities obtained by analyzing the published data from Montes 
et al. (2001) are also accurate. 

(3) As our statistical simulations show, for all six stars from Table 2, the probability 
of their penetration into the Oort cloud region is essentially zero. 

3.2 Stars with Indirect Distance Estimates 

We do not set the objective to survey all stars with spectrophotometric distance estimates. 
Note only two interesting stars, SSSPM J1549-3544 and SDSS J1416±1348, which we have 
been able to reveal using data from Scholz et al. (2004) and Schmidt et al. (2010a; 2010b). 

(1) According to Scholz et al. (2004), the kinematic data for SSSPM J1549-3544 are: 
a = 15 A 48 m 40 s .23, 5 = -35°44'25".4, /i a cos5 = -591 ± 8 mas yr" 1 , fi 5 = -538 ± 5 mas 
yr _1 , and d spec = 3-4 pc with an error of 1 pc. This star may be the single cool white 
dwarf closest to the Sun (closer than the well-known Van Maanen star, d = 4.3 pc). 
Scholz et al. (2004) obtained the distance d spec = 4 ± 1 pc using 2MASS photometry and 
dspec = 3 ± 1 pc using less accurate SSS (SuperCOSMOS Sky Survey) photometry. Since 
there are no radial velocity data, the space velocities (U, V, W) of this star were estimated 
by Scholz et al. (2004) for three model radial velocities, -50, 0, and 50 km s _1 . 

Following Scholz et al. (2004), we will use d — 4 ± 1 pc (e^/d = 25%) and various 
model radial velocities for this star. Several trajectories are shown in Fig. 2. As we see 
from the figure, the sign of the radial velocity determines whether the star could encounter 
the Sun in the past (positive sign) or in the future (negative sign). We clearly see that 
the encounter with the Sun becomes increasingly close with increasing magnitude of the 
radial velocity. Because of the distance estimation error ±1 pc, formally there is a chance 
of very close encounters. For example, when the plot is shifted vertically downward by 
ml pc, the stellar trajectories highlighted by the heavy lines fall at the boundary of the 
Oort cloud. 

Table 3 gives two results obtained at fixed radial velocities of the star with the addition 
of random errors to the proper motion components and the distance d. We constructed 
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10 000 model orbits for each case and in ^900 cases the star falls into the Oort cloud 
region, d min < 0.485 pc, then P 1 = 900/10000. 

(2) According to Schmidt et al. (2010a, 2010b), the kinematic data for SDSS 
J1416+1348 are: a = 14 /l 16 m 24 s .08, 5 = 13°48'26".7, pi a cos5 = 88.0 ± 2.8 mas yr" 1 , 
Us = 139.9 ± 1.3 mas yr -1 , and V r = —42.2 ± 5.1 km s -1 . The most reliable distance 
estimate, d — 8.0 ± 1.6 pc (e d /d = 20%), was derived by Schmidt et al. (2010a) as a 
mean of five photometric and spectroscopic determinations from 2MASS infrared data 
and SDSS optical data. 

(3) Our statistical simulations show that both SSSPM J1549-3544 and SDSS 
J1416+1348 have a nonzero probability of penetrating into the Oort cloud region: 
Pi = 0.09 and P 1 = 0.05, respectively (the last column of Table 3). Both stars are 
of great interest in the problem being solved. 

We may conclude that determining the trigonometric parallaxes and radial velocities 
of these stars is topical. 

(4) Let us now turn to the list of 25 stars from the catalog by Gliese and Jahreifi (1991) 
that were revealed by Mullari and Orlov (1996) as candidates for a close encounter with 
the Solar system. Since 16 of them are Hipparcos stars, the results of the analysis of their 
trajectories are presented in Garcia-Sanchez et al. (2001) and Bobylev (2010). 

Three stars are located within d < 10 pc; these are GJ 905, GJ 473, and GJ 3166. Since 
the remaining six stars are farther than 10 pc, they were not included in the RECONS list. 
We have already discussed the encounter parameters of GJ 905 in the previous section. 

The star GJ 3166 (designated as No 456 NN in Mullari and Orlov, 1996) is of great 
interest, because, according to the estimate by Mullari and Orlov (1996), it can encounter 
the Sun to a record distance d m i n = 0.16 pc at time t m i n = 1600 thousand years. How 
ever, as it turned, there is no information about its proper motion in the catalog by Gliese 
and Jahreifi(1991). Therefore, its space velocities U,V,W were calculated by assuming 
the tangential velocity to be zero, V t = 0. 

Taking d = 20.8 pc (photometric distance), /x t = mas yr~ x , and V r = —12 km s _1 for 
GJ 3166, as in the catalog by Gliese and Jahreifi(1991), we find the encounter parameters 
dmin = 0.18 pc and t min = 1708 thousand years, which confirm the result by Mullari and 
Orlov. For the same distance and radial velocity but taking the proper motion components 
fi a cos5 = — 103.3 ± 8. 1 mas yr -1 and {i$ = — 67.7±9.2 mas yr -1 from the UCAC3 catalog 
(Zacharias et al., 2009), we find completely different encounter parameters, d m i n = 15 
pc and tmin = 590 thousand years, which are already less interesting in our problem. 
Note that the absolute proper motions of this star are also available in the XPM catalog 
(Fedorov et al., 2009): fi a cos5 = —87.2 mas yr -1 and fxs = —68.1 mas yr^ 1 . These were 
determined by comparing the 2MASS and Palomar Observatory Sky Survey positions 
referenced to galaxies with a mean error of about 6 mas yr -1 in each coordinate (Bobylev 
et al., 2010). 

Finally, according to the estimate by Mullari and Orlov (1996), GJ 473 can encounter 
the Sun very closely, d m in = 0.29 pc, at time t m i n = 7.5 thousand years. These parameters 
were obtained using the radial velocity V r = —553.7 km s _1 (Gliese and Jahreifi, 1991). 
GJ 473 (LHS 333=FL Vir=Wolf 424 AB) is a close binary system with a known orbit 
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(Torres et al., 1999). The measurements by Tinney and Reid (1998) performed with a 
high-resolution spectrometer yield the system's heliocentric radial velocity V r = 0.9 ±1.7 
km s _1 . Note that previous measurements also gave a low radial velocity for this system: 
V r = —5 ± 5 km s _1 (GCRV, Wilson 1953). The radial velocity in the compilation by 
Gliese and Jahreifi(1991) is probably erroneous. Since the encounter parameters of GJ 
473 calculated with its new radial velocity V r = 0.9 ±1.7 km s _1 , d min = 6 ± 5 pc and 
tmin = — 3 ± 6 thousand years, are no longer the close encounter parameters, this star was 
not included in Table 2. 

4 CONCLUSIONS 

Based on currently available kinematic data, we searched for stars that either encountered 
or will encounter the solar neighborhood within less than 3 pc. We considered stars outside 
the Hipparcos list. For each of them, there is an estimate of its trigonometric parallax with 
a relative error e n /ir < 2%, radial velocity, and proper motion components. We found 
six such stars that are an important supplement to the list of Hipparcos stars closely 
encountering the Solar system (Garcia-Sanchez et al., 2001; Bobylev, 2001). 

For the first time, two single stars, GJ 3379 (G 099-049) and GJ 3323 (LHS 1723), 
have been identified as candidates for a close encounter with the solar orbit. 

We confirmed the remarkable result by Miillari and Orlov (1996) that the star GL 
905 could encounter the Sun fairly closely: d min = 0.93 ± 0.01 pc and t min = 37.1 ± 0.2 
thousand years. 

Two unique stars are located in the immediate solar neighborhood (d < 10 pc)— 
the white dwarf SSSPM J1549-3544 and the L dwarf SDSS J1416±1348. Our statistical 
simulations showed that both of them have a nonzero probability of penetrating into the 
Oort cloud region: Pi = 0.09 and Pi = 0.05, respectively. Determining the trigonometric 
parallaxes and radial velocities for these stars is topical. This task can be accomplished 
using both ground-based and space (e.g., GAIA) observations. 

Based on new data, we showed that the stars GJ 473 and GJ 3166 are not suitable 
candidates for a very close encounter with the Solar system as was presumed previously. 
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Fig. 1. Trajectories of the stars relative to the Sun: (1) GJ 866, (2) GJ 406, (3) GJ 
65, (4) GJ 3323; for GJ 905 and GJ 3379, we give the model trajectories computed by 
taking into account the random errors in the observational data (1000 realizations). The 
trajectories hatch the 3a confidence regions, the Oort cloud region is shaded, the dotted 
line indicates the trajectory of Proxima Cen. 
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Fig. 2. Model trajectories of the star SSSPM J1549-3544 relative to the Sun: d = 4 
pc at V r = +50 km s" 1 (1), V r = +25 km s" 1 (2), V r = -50 km s" 1 (3), and V r = -25 
km s^ 1 (4). 
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Table 1: Stars and multiple systems with trigonometric parallaxes 



Star 


a J2000; 
<5j2000 


fi a , mas yr 1 
fis, mas yr -1 


n, mas 


V r , km s 1 


Ref* 


GJ 65 AB 


l ft 39 m 01 s .54 
-17°57'01".8 


3296 ± 5 
563 ±5 


373.7 ±2.7 


29.0 ±4.6 


e 


GJ 299 


8Hl m 57 s .58 
8°46'22".l 


1099 ± 8 
-5123 ±8 


146.3 ±3.1 


-35 ±5 


e 


GJ 388 AB 


10 h 19 m 36 s .28 
19°52 / 12".l 


-502 ± 8 
-43 ±8 


204.6 ± 2.8 


11.6 ±0.3 


e 


GJ 406 


10 h 56 m 28 s .86 
7°00'52"8 


-3842 ± 8 
-2727 ± 8 


419.1 ±2.1 


19.5 ±0.4 


a 


GJ 473 AB 


12 ft 33 m 17 s .38 
8°46'22".l 


-1730 ±8 
203 ±8 


227.9 ±4.6 


0.9 ± 1.7 


d 


GJ 866 ABC 


22 h 38 m 33 s .73 
-15°17'57".3 


2314 ±5 
2295 ± 5 


289.5 ±4.4 


-50 ±1 


c 


GJ 905 


23 ft 41 m 54 s .99 
44°10'40".8 


111 ±8 
-1584 ±8 


316.0 ± 1.1 


-78.0 ±0.4 


a 


GJ 1111 


8 h 29 m 49 s .35 
26°46'33"7 


-1112 ±5 
-611 ±4 


275.8 ±3.0 


9.0 ±0.5 


b 


GJ 3323 


5 h 01 m 57 s .47 
-6°56'45".9 


-550 ± 5 
-533 ± 5 


187.9 ± 1.3 


42.0 ±0.1 


b 


GJ 3379 


6 h 00 m 03 s .50 
2°42'23".67 


311 ±3 
-42 ±3 


190.9 ± 1.9 


30.0 ±0.1 


b 



Note. The radial velocities were taken from the following papers: (a) Nidever (2002), (b) Montes 
et al.( 2001), (c) Delfosse et al. (1999), (d) Tinney and Reid (1998), (e) Barbier-Brossat and 
Figon (2000); the parallaxes are given according to the RECONS list of 100 nearest stars. 
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Table 2: Candidates for a close encounter with the Solar system 



Star 


SP 


M/Mq 


^min i P C 


t min , thousand yr 


GJ 905 


M5.5V 


0.12 


0.93 ±0.01 


37.1 ±0.2 


GJ 3379 


M3.5V 


0.19 


1.32 ±0.03 


-163 ±3 


GJ 65 AB 


M5.5V/M6.0V 


0.11/0.10 


2.21 ±0.12 


-29 ±2 


GJ 406 


M6.0V 


0.09 


2.24 ±0.02 


-15 ± 1 


GJ 3323 


M4.5V 


0.15 


2.25 ±0.05 


-104 ±3 


GJ 866 ABC 


M5.0V/-/- 


0.11/0.11/0.10 


2.52 ±0.07 


32.3 ±0.3 



Table 3: Stars with spectroscopic or photometric distance estimates 



Star 


SP 


M/Mq 


drain i 


t-mini 


V r , 


Pi 








pc 


thousand yr 


km s _1 




SSSPM J1549-3544 


>DC11 


« 0.5 


1.21 ±0.58 


-72 ±6 


±50 


0.09 


SSSPM J1549-3544 






1.21 ±0.58 


±72 ±6 


-50 


0.09 


SDSS J1416±1348 


L5V 


« 0.08 


1.24 ±0.65 


186 ± 44 


-42.2 ±5.1 


0.05 
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